Phenomenology of ageing in the Kardar-Parisi-Zhang equation 
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We study ageing during surface growth processes described by the one-dimensional Kardar-Parisi- 
Zhang equation. Starting from a flat initial state, the systems undergo simple ageing in both 
correlators and linear responses and its dynamical scaling is characterised by the ageing exponents 
a = —1/3, b — —2/3, Ac = Xr = 1 and z = 3/2. The form of the autoresponse scaling function is 
well described by the recently constructed logarithmic extension of local scale-invariance. 
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The study of the motion of interfaces continues as 
a widely fascinating topic of statistical physics. One 
particularly intensively studied case is non-equilibrium 
growth processes, which are governed by local rules. Of 
these, the model equation proposed by Kardar, Parisi 
and Zhang (KPZ) [1[ continues to play a paradigmatic 
role in the investigation of the dynamical scaling of such 
interfaces, with an astounding range of applications, in- 
cluding Burgers turbulence, directed polymers in a ran- 
dom medium, glasses and vortex lines, domain walls and 
biophysics, see 0-0 for reviews. Remarkably, in ID the 
height distribution can be shown to converge for large 
times towards the gaussian Tracy- Widom distribution 
[lol 11 1- A particularly clean experimental realization 
of this universality class has been found recently in the 
growing interfaces of turbulent liquid crystals [12| . 

New insight in the non-equilibrium properties of many- 
body systems comes from an analysis of the ageing prop- 
erties, which is realised if the system is rapidly brought 
out of equilibrium by a change of one of its state variables 
13L 14 1 . By definition, an ageing system (i) undergoes 
a slow, non-exponential relaxation towards its stationary 
state(s), (ii) does not satisfy time-translation-invariancc 
and (iii) shows dynamical scaling. Studies of ageing re- 
quire the analysis of both correlators C and responses R 
to be complete and also go beyond the study of dynamics 
in analysing at least two-time observables. Let s denote 
the waiting time and t > s the observation time. For sim- 
ple ageing, one expects in the ageing regime s T m j cro 
and t — s ^> T m j cro , where r m i cro is a microscopic time 
scale, a single relevant length scale L(t) ~ t x / z such that 
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where j is the external field conjugate to </>. This de- 
fines the ageing exponents a, b and the scaling functions, 



from whose asymptotic behaviour fc,B.{y) ~ y~ Xc < R / z as 
y — > 00 one has the autocorrelation and autoresponse 
exponents Ac. r where z is the dynamical exponent. In 
the context of Janssen-de Dominicis theory, cj>(t) is the 
response field conjugate to the order-parameter (j>(t). 

For example, simple ageing is found in non-disordered, 
unfrustrated magnets, quenched from an initial disor- 
dered state to a temperature T < T c at or below its 
critical temperature T c (see [lij and refs. therein) or 
else in microscopically irreversible systems with a non- 
equilibrium stationary state flB - fli! ]. Generically, one 
finds Ac = Xr, but the values of a, b depend more sen- 
sitively on the kind of ageing investigated (for reversible 
systems on the type of quench and for irreversible ones 
on the specific type of dynamics). 

Here, we shall study what kind of ageing phenomena 
can arise in the growth of interfaces. A typical system 
is formulated in terms of a height variable h = hi(t) = 
h{t, Yi), defined over a substrate in d dimensions. A local, 
microscopic rule indicates how single particles are added 
to the surface. One of the main quantities studied is the 
surface roughness 



u 2 (t;L) = ± i J2((h i (t)-h(t)) 2 



(2) 



on a lattice with L d sites and average hei ght h(t) 
L~ d J2i hi(t). It obeys Family- Vicsek scaling [lj| 
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where j3 is the growth exponent and £ = j3z is the rough- 
ness exponent. For an infinite system, the width grows 
for large times as w 2 {t; 00) ~ t 2f) . 

The generic universality class for growth phenomena is 
given by the KPZ equation jl[ 
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TABLE I: Some dynamical, ageing and growth exponents of 
several universality classes in d — 1 dimension. 



where 77 (i, r) is a white noise with zero mean and vari- 
ance (r)(t,r)r)(t',r')) = 2uT6(t - t')6(r - r') and n,v,T 
are material-dependent constants. For comparison, we 
introduce two more universality classes of surface growth: 
elimination of the non-linear term in (0| by setting fi = 
gives the Edwards- Wilkinson (EW) universality class 
[20(. The Mullins- Herrings (MH) universality class is 
given by d t h = -vd*h + r\ (21]. For both EW and MH 
classes, the ageing scaling forms ([T]) for C and R have 
been explicitly confirmed [22|. Values of some growth 
and ageing exponents in II? are listed in table HI For the 
ID KPZ class, the exponents z,f3 are exactly known [l|, 
whereas the relation 6=— 2£/z = — 2/3 follows from dy- 
namical scaling 123L | 24J . Also, evidence for a growing 
length L(t) ~ ^/^pHll and estimates of A c [H, [24[ 
have been reported [27|. However, no systematic test of 
the ageing scaling has been reported for the space-time 
correlation function, and no information exists at all for 
the response function R. These will be provided now. 

Our numerical simulations in the ID KPZ class ei- 
ther use the discretised KPZ equation (j4]) in the strong 
coupling limit [28| (we checked that our results do not 
depend on the chosen discretisation scheme) or else 
the Kim-Kosterlitz (KK) model (2^|. This model uses 
a height variable hi(t) £ Z attached to the sites of 
a chain with L sites and subject to the constraints 
\hi(t) — hi±i(t)\ = 0,1, at all sites i. From a fiat initial 
condition, that is /i;(0) = 0, the dynamics of the model 
is as follows: at each time step, select randomly a site 
i and deposit a particle with probability p or else elimi- 
nate a particle with probability 1 — p. L such deposition 
attempts make up a Monte Carlo step. It is well-known 
that this model is in the KPZ universality class. The 
choice of the value of p is a practical matter. In order 
to avoid meta-stable states, we have chosen p = 0.98. 
In simulations, we have taken L = 2 17 and all the data 
have been averaged over 10 4 samples. For the discretised 
KPZ equation we considered systems of size L = 10 4 and 
averaged over typically 10 5 samples. 

In studying the ageing behaviour, we shall consider the 
two-time spatio-temporal correlator 

C{t,s;v) = ((h(t,r + r )-(h(t)))(h(s,r )-(h(s)))) 
= (h(t,r + r )h(s,r )) - (h(t)) (h(s)) 

(5) 
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FIG. 1: Scaling of the two-time autocorrelator C(t,s) from 
the KK model and several values of s, and the discretised 
KPZ equation, with s — 400, rescaled by a factor 2.79. 



L — > 00 limit and where the definition of the exponents 
is analogous to the usual one for simple ageing. The 
autocorrelation exponent can be found from fc(y) = 
Fc(y,0) ~ y~ x °l z as y -> 00. We also have b = -2/3, 
since the width w 2 (t\ 00) = C(t, t; 0) = t~ h Fc{l, 0). This 
is justified since the initial conditions in the ID KPZ do 
not generate new, independent renormalisations [30| . 

In figure [TJ we show data for the autocorrelator 
C(t, s) = C(t, s; 0) obtained from the KK model. A clear 
data collapse is seen and for large values of the scaling 
variable y = t/s, an effective power-law behaviour with 
an exponent Xc/z « | is found. The data are fully com- 
patible with a numerical solution of the KPZ equation 
and directly test simple ageing ([T]) in the ID KPZ class. 
All this, completely analogous to the EW and MH classes, 
confirms and strengthens earlier conclusions 23-2(1 3pj . 

In order to define a response, we appea l to the proce- 



dures used in irreversible systems [13 1121 yj, 1311 where 
the external field is related to a local change of rates. 
In the KK model, we consider a space-dependent depo- 
sition rate Pi = pa + a^e/2 with a, = ±1 and e = 0.005 a 
small parameter. Then consider, with the same stochas- 
tic noise 77, two realisations: system A evolves, up to 
the waiting time s, with the site-dependent deposition 
rate pi and afterwards, with the uniform deposition rate 
Pq. System B evolves always with the uniform deposition 
rate pi = pq. Of course, the evaporation rate g, = 1 —pi- 
Then, the time-integrated response function is 



X(t,s;r) = / duR(t,u;r) 
Jo 
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along with the extended Family- Vicsek scaling in the 



together with the expected scaling. The time-integrated 
autoresponse x(i, s) = x(t, s; 0) plays the same role as 
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FIG. 2: Scaling of the two-time integrated autoresponse 
x(t, s) from the KK model for several values of s, and from 
the KPZ equation with s = 400, rescaled by a factor 7.25. 



the thermoremanent integrated response of magnetic sys- 
tems Q. The autoresponse exponent is read off from 



oo. For the discre- 



f x (y) = F x (y,0) ~ y-^' z for y 
tised KPZ equation we realize the perturbation by adding 
a spatially random force, of strength ±/q = ±0.3, up to 
the waiting time s. 

In figurc[2] data for the integrated autoresponse x(t, s) 
coming from the KK model are shown. An excellent col- 
lapse is found for a = — |. The effective power law, for 
y = t/s large, reproduces well the expected Xr/z w |. 
Indeed, from the exact fluctuation-dissipation theorem 
TR{t,s;r) = -d*C(t,s;r), valid in the ID KPZ univer- 
sality class (because of time- reversal- invariance) (32-35[. 
we obtain the predictions l + a = b + 2/z and Ac = Xr, 
in agreement with our data. The data are essentially 
identical to those obtained from the KPZ equation, in 
agreement with universality. For the first time, the age- 
ing form ([T|) of the linear response is confirmed in a non- 
linear growth model. In contrast with the EW and MH 
classes, a and b are different, a feature commonly seen in 
irreversible systems [3, . 

Next, in figure [3l we illustrate the space-dependent 
scaling of both the correlator and the integrated response. 
For several values of the scaling variable y — t/s, the de- 
pendence on the second argument in the scaling forms 
(|5l6p is illustrated. An excellent data collapse is found, 
which further confirms the conclusions already drawn 
from the autocorrelator and the autoresponse and also 
confirms the exactly known dynamical exponent z = § 
in the ID KPZ universality class. The shape of the scal- 
ing functions changes notably when y is varied. 

We now turn to an analysis of the form of the autore- 
sponse scaling function f x {y). For ageing simple magnets 
(i.e. non-disordered and unfrustrated) , it has been pro- 
posed to generalise dynamical scaling to a larger set of 
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FIG. 3: Space-dependent scaling of the two-time correlator 
C(t, s) (upper row) and the integrated response x{t, s) (lower 
row), for several values of s and the scaling variable y = t/s. 



local scale-transformations |37j , which includes the trans- 
formation 1 1-> t/(l + jt). This hypothesis of local scale- 
invariance (LSI) indeed reproduces precisely the univer- 
sal shapes of responses and correlators in a large variety 
of models, as reviewed in detail in An alog ous evi- 

dence exists in some irreversible models [3, HI 17. 36j|. 
Similarly, the responses and correlators in the EW and 
MH classes, with the the local height variable h(t, r)—h(t) 
acting as a quasi-primary scaling operator, are described 



by LSI [22 



|. Is LSI also realised in the ID KPZ class, 
with the local height as a quasi-primary operator ? 

We concentrate on the autoresponse function and shall 
restrict attention to the transformations in time. The 
transformation 8(f) = eX n 4> of the quasi- primary oper- 
ators under local scale-transformations is given by the 
infinitesimal generators X n , which read [39j ] 
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and satisfy the commutator [X n ,X„ 



, n > (7) 

(n - m)X n+m . 
We merely look at the finite-dimensional sub-algebra 
spanned by the dilatations Xq and the special transfor- 
mations X\. Since time-translations (generated by A_i) 
are absent, we have two distinct scaling dimensions x 
and £, which together give the shape of the autoresponse 
function, see below. Now, consider a possible extension 
to so-called logarithmic form, where a primary operator 



is replaced by a doublet 







In analogy with log- 



4ol 411, this extension is 



arithmic conformal invariance 
formally carried out by replacing the scaling dimensions 
x, £ by matrices (restricted to the 2x2 case) [42[ 
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where the first scaling dimension is immediately taken in 
a Jordan form. Consistency with the commutators then 
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leads to £" = [42[ . Recalling ([T]) , consider the following 
quasi-primary two-point functions, with y = t/s 
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where .F(y) = y 
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plicitly known scaling functions |42| . In contrast to loga- 
rithmic conformal invariance, logarithmic corrections to 
scaling are absent if x' = x' = and there are no loga- 
rithmic factors for y — > oo if furthermore £' = 0. If we 

take R(t,s) = (ifj(t)^(s)\ = s _1-0 /.R(*/ s )i we find 



fR(y) = y- XR/z (i-y- 1 ) 
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with the exponents 1 + a = + a' — a = | (£ + £ 

Xr/z = x + £ and the normalisation constants ho, go, fo- 
The integrated autoresponse x(t, s) = s~ a f x (t/s) is 
found from (|10j) by using the specific value Xr/z — a = 1 
which holds true for the ID KPZ. We find 

f x (y) = y +1/3 {A [i-(i-y- i y a '} (ii) 

+ (1 - y- 1 ) a ' [Ai In (1 - y- 1 ) + A 2 In 2 (l - y" 1 )] } 

where ^4o,i.2 are normalisations related to fo,go, ho. The 
non-logarithmic case is recovered if A\ = ^2=0. Indeed, 
for y 1, one has / x (y) ~ y~ 2 ^ 3 , as it should be. 

In figure HI which gives a more fine apprecia- 
tion of the shape of f x (y) than figure [2j we com- 
pare data for the reduced scaling function f re d(y) = 
f x (y)y^ 1 ^ 3 [l - (1 - y -1 ) 1 / 3 ] with the predicted form 
(fTTj). Data with s < 10 3 are not yet fully in the scaling 
regime. If one tries to fit the data with a non-logarithmic 
LSI (then R = ((fxj)) or (ip<f>)) one obtains an agreement 
with the data, with a numerical precision of about 5%. 
An attempt to fit only with the first-order logarithmic 
terms (then R = (0^)) with — assumed, gives back 
the same result, see table HH Only if one uses the full 
structure of logarithmic LSI, an excellent representation 
of the data is found, to an accuracy better than 0.1% 
over the range of data available. In the inset the ra- 
tio x(ij s )/Xl 2 lsi(^> s ) i s shown and we see that at least 
down to t/s fa 1.03, the data collapse indicating dynam- 
ical scaling holds true, within the accuracy limits set by 
the stochastic noise, within « 0.5%. For the largest wait- 
ing time s = 4000, this observation entends over the en- 
tire range of values of t/s considered. This indicates that 
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FIG. 4: Comparison of the reduced scaling function /red(y) = 

fx(y)y- 1/3 (i - (i - i/y) 173 )- 1 of x (t,«) = s^f x (t/s) with 

logarithmic local scale-invariance. Non-logarithmic LSI gives 
the dash-dotted curve labelled LSI and full logarithmic LSI 
(|lip gives the dashed curve labelled L 2 LSI. The inset shows 
the ratio x(t, s )/Xl 2 lS!^' s ) over a S ams t */ s — 1- 
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R 


a 1 A Ai A 2 


(#) - LSI 
{<$) - L X LSI 
tyi)) - L 2 LSI 


-0.500 0.662 
-0.500 0.663 -6 • 10^ 4 
-0.8206 0.7187 0.2424 -0.09087 



TABLE II: Fitted parameters j4q,i,2 and a' used in figure [4] 



the local height h and its response field h of the ID KPZ 
equation could be tentatively identified with the logarith- 
mic quasi-primary operators ip, tp, which slightly gener- 
alises the findings for the EW and MH classes, which obey 
non- logarithmic LSI. It is an open question whether the 
approach outlined here just generates the first two terms 
of an infinite logarithmic series in R(t, s). 

A systematic analysis of the invariance properties of 
the dynamical functionals studied for instance in [30l 33 . 
13411 . or the alternate form derived in 431, would be of 



interest, following the lines of study for the analysis of 
dynamical symmetries in phase-ordering kinetics, see fldj 
and refs. therein. 

Summarising, we tested the full scaling behaviour of 
simple ageing, both for correlators and responses, of sys- 
tems in the ID KPZ universality class. This is the first 
example of a growth process described by a non-linear 
equation which is shown to satisfy simple ageing scaling 



5 



for space- and time-dependent quantities. It is non-trivial 
that the values of the growth and dynamical exponents, 
previously known from the study of the stationary state, 
are confirmed far from stationarity. It would be inter- 
esting to measure it also experimentally. Performing a 
numerical experiment, we find the form of the autore- 
sponse scaling function to be very well described by the 
recently constructed logarithmic extension of local scale- 
invariance, with a natural identification of the leading 
quasi-primary operators. In view of important recent 
progress in the exact solution of the ID KPZ equation, 
see p], EH [H| , one may expect that the question of a loga- 
rithmic dynamical scaling can be adressed and its further 
consequences explored. 
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